Statistics Testing

Anova test

sq_df = read_csv('squirrel_tidy.csv') %>% 
  janitor::clean_names()

long_model_final = lm(long ~ shift + age + activity + reaction + sounds, data = sq_df)
anova(long_model_final)

## Analysis of Variance Table
## 
## Response: long
##             Df   Sum Sq    Mean Sq F value    Pr(>F)    
## shift        1 0.000250 0.00025024  4.2452  0.039446 *  
## age          1 0.000446 0.00044617  7.5690  0.005973 ** 
## activity     1 0.000301 0.00030085  5.1038  0.023944 *  
## reaction     1 0.000423 0.00042295  7.1751  0.007432 ** 
## sounds       1 0.001139 0.00113910 19.3243 1.141e-05 ***
## Residuals 3017 0.177841 0.00005895                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

lat_model_final = lm(lat ~ primary_fur_color + reaction + sounds, data = sq_df)
anova(lat_model_final)

## Analysis of Variance Table
## 
## Response: lat
##                     Df   Sum Sq   Mean Sq F value    Pr(>F)    
## primary_fur_color    1 0.001256 0.0012564  12.034 0.0005297 ***
## reaction             1 0.001134 0.0011339  10.861 0.0009936 ***
## sounds               1 0.002101 0.0021015  20.128 7.512e-06 ***
## Residuals         3019 0.315195 0.0001044                      
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1